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ABSTRACT. We introduce a new and highly tractable structural model for spot and derivative prices 
in electricity markets. Using a stochastic model of the bid stack, we translate the demand for power 
and the prices of generating fuels into electricity spot prices. The stack structure allows for a range of 
generator efficiencies per fuel type and for the possibility of future changes in the merit order of the 
fuels. The derived spot price process captures important stylized facts of historical electricity prices, 
including both spikes and the complex dependence upon its underlying supply and demand drivers. 
Furthermore, under mild and commonly used assumptions on the distributions of the input factors, we 
obtain closed-form formulae for electricity forward contracts and for spark and dark spread options. As 
merit order dynamics and fuel forward prices are embedded into the model, we capture a much richer 
and more realistic dependence structure than can be achieved by classical reduced-form models. We 
illustrate these advantages by comparing with Margrabe's formula and a simple cointegration model, 
and highlight important implications for the valuation of power plants. 

Electricity markets and structural model and forward prices and spread options and power plant 
valuation 
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1. Introduction 

Since the onset of electricity market deregulation in the 1990s, the modeling of prices in these 
markets has become an important topic of research both in academia and industry. Energy companies 
own large portfolios of generation units and require sophisticated models of price dynamics in order 
to manage risk. Asset valuation is also of utmost importance in capital intensive industries and 
real option theory is typically used to associate the management of a plant to a string of spread 
options, spanning many years or decades. One of the main thrusts of this paper is to provide new 
and versatile tools for these valuations which efficiently capture the complex dependencies upon 
demand and production fuel prices. 

In electricity price modeling, important challenges include prominent seasonalities and mean- 
reversion at various time scales, sudden spikes and the strong link between the prices of electricity 
and other energy commodities (see Figures lb and 5 for sample daily historical spot and forward 
prices from the PJM market) - features which mostly stem from the non-storability of electricity 
and the resulting matching of supply and demand at all times. While every model should attempt 
to capture these properties as well as possible, at the same time there is a need for fast and efficient 
methods to value power plants and other derivatives on the spot price. To achieve the latter goal 
of efficiency much literature has ignored or oversimplified the former goal of modeling structural 
relationships. In this paper, we propose a model that realistically captures the dependency of power 
prices on their primary drivers; yet we obtain closed-form expressions for spot, forward and option 
prices. 
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The existing literature on electricity price modeling can be approximately divided into three cat- 
egories. At one end of the spectrum are so called full production cost models. These rely on knowl- 
edge of all generation units, their corresponding operational constraints and network transmission 
constraints. Prices are then typically solved for by complex optimization routines (cf. lfT9l ). Al- 
though this type of model may provide market insights and forecasts in the short term, it is — due to 
its complexity — unsuited to handling uncertainty, and hence to derivative pricing or the valuation 
of physical assets. Other related approaches which share this weakness include models of strategic 
bidding (cf. [22J) and other equilibrium approaches (cf. (6)). At the other end of the spectrum are 
reduced form models. These are characterised by an exogenous specification of electricity prices, 
with either the forward curve (cf. [12| and |5||) or the spot price (cf. [16, 20 4|) representing the 
starting point for the model. Reduced form models typically either ignore fuel prices or introduce 
them as exogenous correlated processes; hence they are not successful at capturing the important 
afore mentioned dependence structure between fuels and electricity. Further, spikes are usually only 
obtained through the inclusion of jump processes or regime switches, which provide little insight 
into the causes that underly these sudden price swings. 

In between these two extremes is the structural approach to electricity price modeling, which 
stems from the seminal work of Barlow Q. We use the adjective structural to describe models, 
which — to varying degrees of detail and complexity — explicitly approximate the supply curve in 
electricity markets (commonly known as the bid stack due to the price-setting auction). The market 
price is then obtained under the equilibrium assumption that demand and supply must match. In 
Barlow's work the bid stack is simply a fixed parametric function, which is evaluated at a random 
demand level. Later works have refined the modeling of the bid curve and taken into account its 
dependency on the available capacity (cf. |E1|7][T0)), as well as fuel prices (cf. llZ7lfT4l fTll2l') and 
the cost of carbon emissions (cf. ll23l[T3l ). The raison d'etre of all structural models is very clear. 
If the bid curve is chosen appropriately, then observed stylized facts of historic data can be well 
matched. Moreover, because price formation is explained using fundamental variables and costs of 
production, these models offer insight into the causal relationships in the market; for example, prices 
in peak hours are most closely correlated with natural gas prices in markets with many gas 'peaker' 
plants; similarly, price spikes are typically observed to coincide with states of very high demand 
or low capacity. As a direct consequence, this class of models also performs best at capturing the 
varied dependencies between electricity, fuel prices, demand and capacity. 

The model we propose falls into the category of structural models. Our work breaks from the cur- 
rent status quo by providing closed-form formulae for the prices of a number of derivative products 
in a market driven by two underlying fuels and featuring a continuum of efficiencies (heat rates). In 
the considered multi-fuel setting, our model of the bid stack allows the merit order to be dynamic: 
each fuel can become the marginal fuel and hence set the market price of electricity. Alternatively, 
several fuels can set the price jointly. Despite this complexity, under only mild assumptions on the 
distribution, under the pricing measure, of the terminal value of the processes representing electric- 
ity demand and fuels, we obtain explicit formulae for spot prices, forwards and spread options, as 
needed for power plant valuation. Moreover, our formulae capture very clearly and conveniently 
the dependency of electricity derivatives upon the prices of forward contracts written on the fuels 
that are used in the production process. This allows the model to easily 'see' additional information 
contained in the fuel forward curves, such as states of contango or backwardation — another feature, 
which distinguishes it from other approaches. 

The parametrization of the bid stack we propose combines an exponential dependency on de- 
mand, suggested several times in the literature (cf. l30l fTTl l24l ). with the need for a heat rate 
function multiplicative in the fuel price, as stressed by Pirrong and Jermakyan [27|. Eydeland and 
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Geman [18| propose a similar structure for forward prices and note that Black-Scholes like deriv- 
ative prices are available if the power price is log-normal. However, this requires the assumption 
of a single marginal fuel type and ignores capacity limits. Coulon and Howison lfl4ll construct the 
stack by approximating the distribution of the clusters of bids from each technology, but their ap- 
proach relies heavily on numerical methods when it comes to derivative pricing. In the work of Aid 
et al 1 1 1, the authors simplify the stack construction by allowing only one heat rate (constant heat 
rate function) per fuel type, a significant oversimplification of spot price dynamics for mathematical 
convenience. Aid et al |2J extend this approach to improve spot price dynamics and capture spikes, 
but at the expense of a static merit order, ruling out, among other things, the possibility that coal and 
gas can change positions in the stack in the future. In both cases the results obtained by the authors 
only lead to semi-closed form formulae, which still have to be evaluated numerically. 

The importance of the features incorporated in our model is supported by prominent develop- 
ments observed in recent data. In particular, shale gas discoveries have led to a dramatic drop in US 
natural gas prices in recent years, from a high of over $13 in 2008 to under $3 in January 2012. Such 
a large price swing has rapidly pushed natural gas generators down the merit order, and highlights 
the need to account for uncertainty in future merit order changes, particularly for longer term prob- 
lems like plant valuation. In addition, studying hourly data from 2004 to 2010 on marginal fuels in 
the PJM market (published by Monitoring Analytics), we observe that the electricity price was fully 
set by a single technology (only one marginal fuel) in only 16.1 per cent of the hours. For the year 
2010 alone, the number drops to less than 5 per cent. Substantial overlap of bids from different fuels 
therefore exists, and changes in merit order occur gradually as prices move. We believe that our 
model of the bid stack adheres to many of the true features of the bid stack structure, which leads to 
a reliable reproduction of observed correlations and price dynamics, while retaining mathematical 
tractability. 

2. Structural approach to electricity pricing 

In the following we work on a complete probability space (i2, J^,P). For a fixed time horizon 
T e K+, we define the (« + 1) -dimensional standard Wiener process (W ( °, W r ), G [ r i, where W := 
(W 1 ,. ■ -,W n ). Let^° :— ( J7 ) denote the filtration generated by W° and & w := the filtration 

generated by W. Further, we define the market filtration J£" := V . All relationships between 
random variables are to be understood in the almost surely sense. 

2. 1 . Price Setting in Electricity Markets. We consider a market in which individual firms generate 
electricity. All firms submit day-ahead bids to a central market administrator, whose task it is to 
allocate the production of electricity amongst them. Each firm's bids take the form of price-quantity 
pairs representing an amount of electricity the firm is willing to produce, and the price at which the 
firm is willing to sell iQ An important part is therefore played by the merit order, a rule by which 
cheaper production units are called upon before more expensive ones in the electricity generation 
process. This ultimately guarantees that electricity is supplied at the lowest possible pricen 

Assumption 1. The market administrator arranges bids according to the merit order and hence in 
increasing order of costs of production. 



Alternatively, firms may, in some markets, submit continuous bid curves, which map an amount of electricity to the 
price at which it is offered. For our purposes this distinction will however not be relevant. 

Irrn 

This description is of course a simplification of the market administrator's complicated unit commitment problem, 
typically solved by optimization in order to satisfy various operational constraints of generators, as well as transmission 
constraints. Details vary from market to market and we do not address these issues here, as our goal is to approximate the 
price setting mechanism and capture the key relationships needed for derivative pricing. 
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We refer to the resulting map from the total supply of electricity and the factors that influence 
the bid levels to the price of the marginal unit as the market bid stack and assume that it can be 
represented by a measurable function 

b: [0,|] xR" 9 K,s)4^,s)£l, 

which will be assumed to be strictly increasing in its first variable. Here, £ € K+ represents the 
combined capacity of all generators in the market, henceforth the market capacity (measured in 
MW), and s€K" represents factors of production which drive firm bids (e.g. fuel prices). 

Demand for electricity is assumed to be price-inelastic and given exogenously by an jF,°-adapted 
process (D t ) (measured in MW). As we shall see later, the prices of the factors of production used 
in the electricity generation process will be assumed to be j^^-adapted. So under the objective 
historical measure P, the demand is statistically independent of these prices. This is a reasonable 
assumption as power demand is typically driven predominantly by temperature, which fluctuates at 
a faster time scale and depends more on local or regional conditions than fuel prices. The market 
responds to this demand by supplying an amount £ [0, £,] of electricity. We assume that the market 
is in equilibrium with respect to the supply of and demand for electricity; i.e. 

(1) D f = &, iorte{0,T}. 

This implies that D t e [0, f ] for t e [0, T] and (&) is ^-adapted. 

The market price of electricity (P t ) is now defined as the price at which the last unit that is needed 
to satisfy demand sells its electricity; i.e. using ([TJ, 

(2) P t :=b{D u -), forte [0,r]. 

We emphasize the different roles played by the first variable (i.e. demand) and all subsequent 
variables (i.e. factors driving bid levels) of the bid stack function b. Due to the inelasticity as- 
sumption, the level of demand fully determines the quantity of electricity that is being generated; all 
subsequent variables merely impact the merit order arrangement of the bids. 

Remark 1. The price setting mechanism described above applies directly to day-ahead spot prices 
set by uniform auctions, as in most exchanges today. However, we believe that in a competitive 
market with rational agents, the day-ahead auction price also serves as the key reference point for 
real-time and over-the-counter prices. 

2.2. Mathematical Model of the Bid Stack. From the previous subsection, it is clear that the 
price of electricity in a structural model like the one we are proposing depends critically on the 
construction of the function b. Before we explain how this is done in the current setting, we make 
the following assumption about the formation of firms' bids. 

Assumption 2. Bids are driven by production costs. Furthermore, 

(1) costs depend on fuel prices and firm-specific characteristics only; 

(2) firms' marginal costs are strictly increasing. 

To back up Assumption|2] we briefly consider it in light of historic bid and spot price data. Figure 



la plots the bid stack from the PJM market region in the US on the first day of two consecutive 



months. Firstly, rapidly increasing marginal costs lead to the steep slope of the stack near the market 



capacity (70,000 MW); this feature translates directly into spikes in spot prices (see Figure lbi 
when demand is high. Between the two sample dates in February and March 2003, the prices of gas 
increased rapidly. In the PJM market gas fired plants have historically featured mainly in the second 
half of the bid stack. Therefore, it is precisely the gas price related increase in production costs, 
which explains the increase in bid levels observed beyond about 40,000MW in Figure [Ta] (see also 



{ 4.2 for a discussion of recent merit order changes). 
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Figure 1 . Historical prices and bids from the PJM market in the North East US 



Production costs are typically linked to a particular fuel price (e.g. coal, natural gas, lignite, oil, 
etc.). Furthermore, within each fuel class, the cost of production may vary significantly, for example 
as old generators may have a higher heat rate (lower efficiency) than new units. It is not our aim 
to provide a mathematical model that explains how to aggregate individual bids or capture strategic 
bidding. Instead, we group together generators of the same fuel type and assume the resulting bid 
curve to be exogenously given and to satisfy Assumption [2] From this set of bid curves, the merit 
order rule then determines the construction of the market bid stack. 

Let / = {1, . . . ,«} denote the index set of all the fuels used in the market to generate electricity. 
We assume that their prices are the only factors influencing the bid stacks. With each i £ I we 
associate an J^" r w -adapted fuel price process (S' t ) and we define the fuel bid curve for fuel i to be a 
measurable function 

ft, : [0, x K 3 eR, 
where the argument £, represents the amount of electricity supplied by generators utilizing fuel type 
i, s a possible value of the price S' t , and £ M + the aggregate capacity of all the generators utilizing 
fuel type i. We assume that bj is strictly increasing in its first argument, as required by Assumption 
[2] Further, also for i £ I, let the ^-adapted process (<!;/) represent the amount of electricity supplied 
by generators utilizing fuel type i. It follows that 

Ll'=l, and D, = £#, forte[0,T}. 

In order to simplify the notation below, for i £ /, and for each seMwe denote by b,( ■ ,s)~ l the 
generalized (right continuous) inverse of the function % bi(%,s). Thus 

b,{-,s)-\p) := |'Ainf{§ G (0,|'] : bfas) > p) , for (p,s) ERxR, 

where we use the standard convention inf0 = +°°. Using the notation 

&(*):= 6i(0,j) and 

and writing bj (p,s) =bj(- ,s) (p) to ease the notation, we see that »■ (p,s) — if p G (—00,^(5)), 
bj\p,s)=l i \ipE[b i {s),oo), i mdb;\p 1 s)E[Q,l i ] if p £ (*),*/(*)). For fuel i £l at price S l t = s 
and for electricity prices below bj(s) no capacity from the ith technology will be available. Similarly, 
once all resources from a technology are exhausted, increases in the electricity price will not lead 
to further production units being brought online. So defined, the inverse function bj 1 maps a given 
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price of electricity and the price of fuel i to the amount of electricity supplied by generators relying 
on this fuel type. 

Proposition 1. For a given vector (D t ,S t ), where S; := (Sj ,.. . ,5"), the market price of electricity 
{P t ) is determined by 

(3) P, = rninfe,. (sj) V sup jp e E : ££r 1 (p, S' t ) < D, j , /or r € [0, r] . 
Proof. By the definition of bj l , the function defined by 

5-W,"-,*"):=£SrW)> 

ie/ 

is, when the prices of all the fuels are fixed, a non-decreasing map taking the electricity price to the 
corresponding amount of electricity generated by the market. Similarly to the case of one fixed fuel 
price, for each fixed set of fuel prices, say s := (*,••• ,s n ), we define the bid stack function | 
b(£,s) as the generalized (left continuous) inverse of the function t, ^-s- b~ l {p,s l ■ ■ ,s n ) defined 
above, namely 

:= winbj (s l ) V sup jp e K : gfor 1 (p,. 1 ) < A , for (§,s) € [0, |] x R", 

where we use the convention sup© = — °°. 

The desired result follows from the definition of the market price of electricity in □ 

2.3. Denning a Pricing Measure in the Structural Setting. The results presented in this paper 
do not depend on a particular model for the evolution of the demand for electricity and the prices 
of fuels. In particular, the concrete bid stack model for the electricity spot price introduced in f|3] 
is simply a deterministic function of the exogenously given factors under the real world measure 
P. However, for the pricing of derivatives in f|4]and Sj5]we need to define a pricing measure Q and 
the distribution of the random factors at maturity under this new measure will be important for the 
results that we obtain later. 

For an ^-adapted process d t , where 8 t := (0,°, Of, . . . , 9"), a measure i? ~ P is characterized 
by the Radon-Nikodym derivative 

(4) dP : = eX H-/o e "- dW "-2Jo I "' dU 
where we assume that (6 t ) satisfies the so-called Novikov condition 



E 



T 2 



exp - / |fl„| z dii 



We identify (0 f °) and (0/) with the market prices of risk for demand and for fuel i respectively. 
We choose to avoid the difficulties of estimating the market price of risk (see for example |[T9l for 
several possibilities) and instead make the following assumption. 

Assumption 3. The market chooses a pricing measure Q ~ P, such that 

Q £ {J2 ~ P : all discounted prices of traded assets are local i?-martingales} . 

Note that we make no assumption regarding market completeness here. Because of the non- 
storability condition, certainly electricity cannot be considered a traded asset. Further, there are 
different approaches to modeling fuel prices; they may be treated as traded assets (hence local mar- 
tingales under Q) or — more realistically — assumed to exhibit mean reversion under the pricing 
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measure. Either way, demand is a fundamental factor and the noise (W t °) associated with it means 
that the joined market of fuels and electricity is bound to be incomplete. Note however, that all 
derivative products that we price later in the paper (forward contracts and spread options) are clearly 
traded assets and covered by Assumption [3] 



3. Exponential Bid Stack Model 

Equation |3]l in general cannot be solved explicitly. The reason for this is that any explicit solution 
essentially requires the inversion of the sum of inverses of individual fuel bid curves. 

We now propose a specific form for the individual fuel bid curves, which allows us to obtain a 
closed form solution for the market bid stack b. Here and throughout the rest of the paper, for i E I, 
we define bi to be explicitly given by 

(5) bi(g,s) := sexpfa+mt), for G [0,|'] x E+, 

where kj and m, are constants and m, is strictly positive. Note that bi clearly satisfies ([l} and since it 
is strictly increasing on its domain of definition it also satisfies (|2j. 

3.1. The Case of n Fuels. For observed (D t ,S t ), let us define the sets M.C C /by 
M := {i E I : generators using fuel i are partially used} 
and C := \i G/ : the entire capacity f ' of generators using fuel i is used} . 

A possible procedure for establishing the members of M and C is to order all the values of b { and bi 
and determine the corresponding cumulative amounts of electricity that are supplied at these prices. 
Then find where demand lies in this ordering. 

With the above definition of M and C we arrive at the following corollary to Proposition [T] 

Corollary 1. For bi of exponential form, as defined in (|5j, the market price of electricity is given 
explicitly by the left continuous version of 

(6) P t =(n(S^yx V h+r\D t -^J i j \, forte[0,T\, 



where 




7\L k i n m J ) » 



and C := m 3' 

leM jeM.j^l 

Proof. At any time t G [0, T] the electricity price depends on the composition of the sets M and C; 
i.e. the current set of marginal and fully utilized fuel types. 

For i E M, Sf 1 = bj x , for i G C, bj l = |' and for i El\{MUC}, br 1 = 0. Therefore, inside 
the supremum in (BJ, we replace / with M and take Y^iec^' to th e right hand side. By Proposition 
[l]the electricity price is given by the left continuous inverse of the function Y.ieM bj , which in the 
exponential case under consideration, simplifies to a single log function and yields |6]). □ 

It is clear from equation |6| that the number of possible expressions for the electricity price is 
fully determined by the different configurations the sets M and C can take. In fact, fluctuations in 
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Supply (MW) Supply (MW) 

(a) Fuel bid curves bj. (B) Market bid stack b. 



Figure 2. Example of fuel bid curves and resulting market bid stack for / := 
{1,2,3},?:= | 



demand and fuel prices can lead to 

n / \ n—i 
M\ / \j=0 

distinct cases for Q. Nonetheless, the market bid stack is always a piece-wise exponential function 
of demand (see Figure[2]i with constantly evolving shape as fuel prices move. This captures in a very 
tractable way the influence of fuel prices on the merit order and resulting power price dynamics. 

3.2. The Case of Two Fuels. For the remainder of the paper, we restrict our attention to the case 
of a two-fuel market, consisting of coal and natural gas generators. Our results can in principle 
be extended to the general case of n > 2 fuels. However, the level of complexity of the formulas 
increases rapidly, as evidenced by the number of possible expressions given in Q. We also choose 
to omit the analysis of the one fuel case, which leads to far simpler expressions throughout, but 
cannot lead to merit order changes. From now on, we set / := {c,g} and carry over all notation 
introduced in f|2]and ^3]so far. 

From |7) we know that there are five possible expressions for the electricity spot price. We list 
them in Table [T] Note that fixing D t reduces this list to some subset of three, each of which — 
depending on the state of S, — can set the electricity price. (A similar reduction to three expressions 
occurs by fixing S,.) We exploit this property to write formula (|6]l in a form more amenable to 
calculations, identifying all cases explicitly. To simplify the presentation in the text below, we 
define 

b cg (£,s) := exp (J3 + 7 §) , for (§,s) G [0, §] x R 2 + , 

where a c , a g , 8 and y are defined in Corollary [T] and simplify for two fuels to 

m e m c „ k c m„ + k„m c m c m„ 

a c = s - — , a„ = l-a c = — , 8 = -^ 8 —, y= — 

m c + m g ' m c + m c + m g m c + m g 

Further, we set /_ := argminjf and i + := argmax{| c ,f*}. 

Corollary 2. With I := {c,g}, for t € [0, T], the electricity spot price is given by 

P, = b, ow (A,S,)I[ ,f.-] (D t )+b mid (D,,S,)I ( |,_|, + ] {D,)+b m (D t ,S t )l^ i+ ^ (A), 
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P t , for t £ [0, T] 


Criterion 


Composition of 
M C 


exp (k c + m c D t ) 


b c (D t ,Sf)<b g (Sf) 


to 


{0} 


Sf exp (kg +m g D t ) 


b g (D t ,Sf)<b c (Sf) 


to 


{0} 


S'exp(k c + m c (D t -¥)) 


b c (D t -¥,%)>b g (S?) 


to 


to 


Sfexp(k g + m g (D t -^)) 


b g (D t -fr,Sf)>b e (S°) 


to 


to 


W) ac (5f) as exp(/3+7A) 


otherwise 


{c,g} 


{0} 



Table 1 . Distinct cases for the electricity price (|6]l in the two fuel case 



where, for (£,s) € [0,^]xl}, 

b hw (4 ,s) := b c (^Ol{ M ^)<i g(sS )} +b g ($,s 8 ) l{b g{ ^ sg)<Us c)} 

+ bc ^^\b c {^)>k g {s«),b g {^)>U^)Y 

b mid {^):=b h tt> si+ )\ bi+{ ^)< h _{j-)} + K tt-?^ si+ )\ bi+ {Z-^A)>h^-)} 

+ beg (S , s) ) ^ (s ,_ ) K ^_ v _ s , + ) (s ,_ ) | , 

high (C.s) := ^ (4 - 1 V) I{* c (|-|*,^)>5 g (5*)} +bg (£ " f V*) I {i g (^| t ',.^)>5 c ( J -)} 

Proof. The expressions for Z?i ow , £> m id, frhigh are obtained from |6} by fixing Z) r in the intervals 
(0,4'~], (<^' + :^] respectively and considering the different scenarios forM and C. □ 




Figure 3 . Illustration of the dependence of power spot price on fuel prices and demand 



Figure [3] illustrates the dependence of the power price on the input factors, namely the prices of 
coal and gas and the demand for electricity. For the purpose of these plots we assumed throughout 
that 4 £ > l^ g . The surfaces in Figure 3a represent the functions b\ ov/ , Z? m ;d and frhigh for demand 



fixed in the corresponding intervals. In Figure 3b we aggregate the information contained in Figure 
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3a for a fixed coal price; i.e. we compare the dependency of the electricity price on the gas price 
for each of the three relevant demand levels. In all three cases electricity is non-decreasing in fuel 
price and is constant against S' t if fuel i is not at the margin (i.e., ; ^ M). Furthermore, in each case 
P t is linear in the sole marginal fuel for sufficiently low or high gas price, and non-linear in both 
fuels in the quadrilateral in the middle. This characterises the region of coal-gas overlap, where both 
technologies jointly set the price. Finally, note that labels PI to P5 indicate regions corresponding 
to rows one to five of Table Q] 

3.3. Extension to Capture Spikes and Negative Prices. In this section we suggest a simple exten- 
sion of the bid stack model in order to more accurately capture the spot price density in markets that 
are prone to dramatic price spikes during peak hours or sudden negative prices at off-peak times. 
Importantly, this modification does not impact the availability of closed-form solutions for forwards 
or spread options, which we introduce in Sj4]and Sj5] 

Let (X t ) be a stochastic process adpated to the filtration j?" r , as is D t . Further, we assume that the 
relationship between the two processes satisfies 

{X t < 0} = {to G Q. : A = 0} and {X, > 0} = { A = I } • 

The difference between (A) and (X t ) is that (A) is restricted to take values in [0,<^] only, whereas 
the process (X,) can potentially take values on the entire real line. 

In the event that demand hits zero or <!;, we say that the market is in a negative price regime 
(for which negative prices are possible, but not guaranteed) or a spike regime and we redefine the 
electricity price at these points to be given respectively by 



b„ (x,s) := b(0,s) — exp(— m„x) + 1, for (x,s) e (— °°,0] 



X 



and b s (x,s) :=fc(f ,s) + exp(/n s (x — f)) — 1, for (x,s) € [f ,°°) x R+- 

Under this extension, the power price expression from Corollary [2] is therefore replaced by 

P t :=b„ {X, , S, ) I {0} (A ) + /'low (A , SO I(o,|*] (A) 

+2>mid (A,S,)I (| ,^(A)+^i gh (A,S,)I (|A|) (A)+/' J (^,S r )I| ?} (A), 

where b\ ow , &mid, &high are defined in Corollary [2] 

Note that the constants m„ , m s > determine how volatile prices are in these two regimes. In such 
cases, the price of electricity may now be interpreted as being set by a thin tail of miscellaneous bids, 
which correspond to no particular technology. Therefore, the difference between the electricity price 
implied by the bid stack and that defined by the negative price or spike regime is independent of fuel 
prices. 

Remark 2. It is possible to generate realistic spikes even in the base model (without the inclusion of 
the spike regime), simply by choosing one of the exponential fuel bid curves to be very steep (large 
mi). However, this would come at the expense of realistically capturing changes in the merit order, 
as it artificially stretches the bids associated with that technology. 

Figure [4] displays the electricity price through time as generated by the stack model for three 



different choices of parameters, for the same scenario. In Figure 4b we show a typical price path in 
the case that m c , m g , m s and m„ are very small. This corresponds to a step function bid stack and 
has been suggested by Aid et al 0]. Clearly, the prices do not exhibit enough variation to match 
observed time series. The solid line in Figure |4a| corresponds to more realistic values of m c and m g ; 
the dashed line illustrates the modification of this path due to the choice of larger values for m s and 
m n . Both paths capture the stylized facts of electricity price time series reasonably well. For the 
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FIGURE 4. Simulation of the power price for typical parameters 

purpose of this simulation the prices of coal and gas have been modeled as exponential Ornstein- 
Uhlenbeck (OU) processes, and demand as an OU process with seasonality truncated at zero and | 
(with (X t ) its untruncated version). However, the choice of model for these factors is secondary at 
this stage, as we are emphasizing the consequences of our choice for the bid stack itself. 

4. Forward Contracts 

We now turn to the analysis of forward contracts in our structural framework. For the sake of 
simplicity, we ignore delivery periods and suggest that T be considered as a representative date in a 
typical monthly delivery period (see [4| for more on handling delivery periods). For the purpose of 
the present discussion, a. forward contract with maturity T is defined by the payoff 

Pr-Ff, 

where F t p is the delivery price agreed at the initial date f , and paid by the holder of the long position 
at T. Simple arbitrage arguments (cf. Il26ll ) imply that 

Ff = E Q [P T \^ t }. 

The result of Corollary [2] shows that the payoff of the forward is a function of demand and fuels, so 
that the electricity forward contract becomes a derivative on fuel prices and demand. 

4.1. Closed Form Forward Prices. For the explicit calculation of forward prices the following 
property of Gaussian densities will be useful (see [21 1 for an application of this result to pricing 
compound options). Let <p\ denote the density of the standard univariate Gaussian distribution, 
and <t>i (•) and ^(t'.P) the cumulative distribution functions (cdfs) of the univariate and bivariate 
(correlation p) standard Gaussian distributions respectively. 

Lemma 1. The following relationship holds between <pi, <t>i and 4>2.' 



(Si / exp(/i + c/i.v up, i.v !4>,(/ 2 c/ 2 .v) dv = exp ( + *2 ( fl-gi,-^=B; 2 



for all li,h,qi,q2 £ K and a £ RU {°°}. 
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Proof. In equation ([8]) combine the explicit exponential term with the one contained in <pi and com- 
plete the square. Then, define the change of variable (x,y) — > (z,w) by 



x = z + qi, y = wyl+qj + q2(x-qi). 

The determinant of the Jacobian matrix J associated with this transformation is |/| = ^Jl + q\. 
Performing the change of variable leads to the right hand side of (|8). □ 

For the main result in this section we denote by F t \ i £ I, the delivery price of a forward contract 
on fuel i with maturity T and write F f := (F t c ,F t c ). 

Proposition 2. Given I = {c,g}, if under Q, the random variables \og(S c T ) and log(Sj.) are jointly 
Gaussian with means jX c and jl g , variances <7 2 and (7 2 and correlation p, and if the demand Dj 
at maturity is independent of ' , then for t £ [0,7*], the delivery price of a forward contract on 
electricity is given by: 

(9) Ff = / f low (D, F, ) fa D) dD+ I f mid D,F, fa (D) dD+ I f high (D, ¥,) fa (D) dD, 

where fa denotes the density of the random variable Dj and, for (4,x) £ [0,4] 
/ W (§, X )=^fc i (§,X i )^l(^ 5 0)/ t 7) 



t>2 . 



/£/ 



- b cg {^x)exp (-a c a g a 2 /2) 



0)/a + a,a) 



f mid ^,x) =b i+ (| -f''V+)<I>i (-* <+ (4 -I*-,!*-) lo) +b l+ (§,0)/a) 
+/7 c . ? (^x)exp(-a t .a ? (J 2 /2) [<J>, - f <-,f '-) /a + a,_cj) -*! (4,0)/a + a ; -_c7)] , 



f m = j> (i - |V) <*>! (-*< (I - !') /a) 



it/ 



-b cg {£,x)exp (-a c a g o 2 /2) 



fa /*+«/<*) 



it/ 



where j = I \ {/}, f/ie constants a t , (X g ,fi ,J are as defined in Corollary^ and 



a 1 := a £ 2 -2p<7 c crg + c7 g 2 , 



Ri := + m;|y - ft* - m^i + log (V/ ) 

Proof. By iterated conditioning, for f £ [0, T], the price of the electricity forward F/ is given by 



-l0g(F/)--CT 2 



(10) 



Ff 



[P T \&,]=^ 



'[b(D T ,S T )\^V^ 



The outer expectation can be written as the sum of three integrals corresponding to the cases Dj £ 
[0,4' _ ], Dj £ [£,'-,£,'+] and Dj £ [4' + ,0] respectively. We consider the first case and derive the /i ow 
term. The other cases corresponding to f m ^ and /high are proven similarly. 

From Corollary [5] we know that b — b\ ow for Dj £ [0, \ l ~ } . This expression for Pj is easily written 
in terms of independent standard Gaussian variables Z := {Zi,Z<i) by using the identity 

^2 



log(^) 
log (if.) 



Mr 



P0c0 g 

2 



Zi 

Z 2 
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Defining d : = E^[Dr|«^r]> the inner expectation can now be written in integral form as 

E[b l0W (d,Z)] =I C +I g +I cg , 

where foi 0W (§,Z) := £>i 0W (<i;,S) and the expectation is computed with respect to the law of Z. For 
example, after completing the square in zi, 

exp(Zi +<7lZ2)0lfe) ( I ) l('2 + <72Z2) dz 2 , 

with 

? Cr ,2 u c + k c + m c d + Up 

2 o c {o c -po g ) 
a g (a g -pa c ) 

q\ :=pa c a g , qi -=-S—± — __. 

a c (a c -pa g ) 

Lemma[T]now applies with a = °°. I g and I cg are computed similarly. In all terms, we substitute for 
Ui using the following standard result. For i £ I, 

(ii) f; = [4 1 & t ] = ex P Ui + l°n, for 1 e [o, r]. 

Substituting the resulting expression for the inner expectation into the outer expectation in ( fT0| ) 
yields the first term in the proposition. □ 



Remark 3. The assumption of lognormal fuel prices in Proposition|2]is a very common and natural 
choice for modeling energy (non-power) prices. Geometric Brownian Motion (GBM) with constant 
convenience yield, the classical exponential OU model of Schwartz 1281 . and the two-factor version 
of Schwartz and Smith |29| all satisfy the lognormality assumption. 

The above result does not depend upon any assumption on the distribution of the demand at 
maturity, and as a result, it can easily be computed numerically for any distribution. In markets where 
reasonably reliable load forecasts exist, one may consider demand to be a deterministic function, in 
which case the integrals in |9| are not needed and the forward price becomes explicit. For cases 
when load forecasts are not reliable, we introduce another convenient special case below, where 
demand at maturity has a Gaussian distribution truncated at zero and f . 

To simplify and shorten the notation we introduce the following shorthand notation: 

,y;pj :=&2{xi,y;p) -Q 2 (x2,y;p)- 
Corollary 3. In addition to the assumptions in Proposition^let demand at maturity satisfy 



D T = max (0,min (t,,X T )) , 
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where Xj ~ N '([14,0%) is independent of St under ' 
forward contract is given explicitly by 



Then for t G [0, T], the delivery price of a 




where j = I\ {/}, <5, = (-1) <''='+> and 



af d :=m?aj + CJ z and rj 

Proof. We use Lemma[T]with a < °o, After integrating over demand, each of the terms in /i ow , /my 
and /hi 2 h turns into the difference between two bivariate Gaussian distribution functions. Simplifying 
the resulting terms leads to the result. □ 

Although the expression in Corollary [3] may appear quite involved, each of the terms can be 
readily identified with one of the five cases listed in Table [T] along with four terms (the last line of 
F t p ) corresponding to the endpoints of the stack. Furthermore, it is noteworthy that as compared to 
Corollary [2] the fuel/orwaraf prices now replace the fuel spot prices in the bid stack curves b{, while 
fl d replaces demand. The Gaussian cdfs essentially weight these terms according to the probability 
of the various bid stack permutations. Thus Ff can become asymptotically linear in Ff or F t g when 
the probability of a single fuel being marginal goes to one. Finally, we note that a very similar 
closed-form expression is available for higher moments of Pf and given in Appendix|A] Convenient 
expressions can also be found for covariances with fuels and for the Greeks (sensitivities with respect 
to underlying factors or parameters), but are not included. 



Remark 4. Under the extended model introduced in £ 3.3 to capture spikes, forward prices are given 
by the same expression as in Proposition[3]plus the following simple terms 



exp m s {ll d -%)- 



2^2 



■m s a d -4> 



1 



exp -m n ll d + -m n O l 




m n a d - 



<3d 



-*1 



Remark 5. In most electricity markets, available capacity is often uncertain, due for example to the 
risk of generator outages. Since enters linearly in the exponential function in |6]) like demand, 
an extension to stochastic capacity levels should be feasible, though rather involved. However, if 
capacity shocks are similar for both fuel types, this additional randomness could more easily be 
accounted for by adjusting demand parameters \i d and a d . More generally, we note that these 
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parameters could in practice be chosen to calibrate the model to observed power forward (or option) 
prices, thus using the random variable Dj as a proxy for demand, capacity, and all other non fuel- 
related risk, along with corresponding risk premia. 

4.2. Correlation Between Electricity and Fuel Forwards. In the American PJM market, coal 
and gas are the fuel types most likely to be at the margin, with coal historically below gas in the 
merit order. Therefore PJM provides a suitable case study for analyzing the dependence structure 
suggested by our model. In Figure [5] we observe the historical co-movement of forward (futures) 
prices for PJM electricity (both peak and off-peak), Henry Hub natural gas (scaled up by a factor of 
ten) and Central Appalachian coal. We pick maturities December 2009 and 2011, and plot futures 



prices over the two years just prior to maturity. Figure 5a covers the period 2007-09, characterized 
by a peak during the summer of 2008, when most commodities set new record highs. Gas, coal 
and power all moved fairly similarly during this two-year period, although the correlation between 
power and gas forward prices is most striking. 
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(A) Dec 2009 forward price dynamics. 
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Figure 5 . Comparison of power, gas and coal futures prices for two delivery dates 



Figure [5b]depicts the period 2009-201 1, during which, due primarily to shale gas discoveries, gas 
prices declined steadily through 2010 and 201 1, while coal prices held steady and even increased a 
little. As a result, this period is more revealing, as it corresponds to a time of gradual change in the 
merit order. Our bid stack model implies that the level of power prices should have been impacted 
both by the strengthening coal price and the falling gas price, leading to a relatively flat power price 



trajectory. This is precisely what Figure 5b reveals, with very stable forward power prices during 
2010-2011. The close correlation with gas is still visible, but power prices did not fall as much as 
gas, as they were supported by the price of coal. Finally, we can also see that the spread between 
peak and off-peak forwards for Dec 2011 delivery has narrowed significantly, as we would also 
expect when there is more overlap between coal and gas bids in the stack. This subtle change in 
price dynamics is crucial for many companies exposed to multi-commodity risk, and is one which 
is very difficult to capture in a typical reduced-form approach, or indeed in a stack model without a 
flexible merit order and overlapping fuel types. 



5. Spread Options 



This section deals with the pricing of spread options in the structural setting presented above. We 
are concerned with spread options whose payoff is defined to be the positive part of the difference 



16 



RENE CARMONA, MICHAEL COULON, AND DANIEL SCHWARZ 



between the market spot price of electricity and the cost of the amount of fuel needed by a particular 
power plant to generate one MWh. If coal (gas) is the fuel featured in the payoff then the option is 
known as a dark (spark) spread. Denoting by h C) h g > the heat rate of coal or gas, dark and spark 
spread options with maturity T have payoffs 

(12) {P T -h c S c T ) + and (P T -h g S*.) + , 

respectively. We only consider the dark spread but point out that all results in this section apply to 
spark spreads if one interchanges c and g. Further, since spread options are typically traded to hedge 
physical assets (generating units) the heat rates that feature in the option payoff are usually in line 
with the efficiency of power plants in the market. Based on the range of market heat rates implied 
by our stack model, we requir^] 

(13) exp (k c ) < h c < exp (k c + m c | e ) . 

Then, as usually, the value (V, ) of a dark spread is given by the conditional expectation under the 
pricing measure of the discounted payoff; i.e. 

V, = exp (-r(T - t))E® [ (P T - h c S c T ) + \ & t ] , 

which thanks to Corollary |2]is understood to be a derivative written on demand and fuels. 

Remark 6. While spread option contracts are often written on forwards, we consider spread options 
on spot prices, as required for our goal of power plant valuation. In addition, as we are interested 
in closed-form expressions, we limit our attention to the payoffs with strike zero, corresponding to 
a plant for which fixed operating costs are negligible or relatively small. Including a positive strike, 
one requires approximation techniques to price a spread option explicitly (such as perturbation of 
the strike zero case), analogously to approaches proposed for when both commodities are lognormal 
(cf. El). 



5.1. Closed Form Spread Option Prices. The results derived in this section mirror those in f4.1 
derived for the forward contract. Firstly, conditioning on demand, we obtain an explicit formula for 
the price of the spread. Secondly, we extend this result to give a closed form formula in the case of 
truncated Gaussian demand. 

We keep our earlier notation for the dominant and subordinate technology i + and i_, and define 



, h _ \ogh c -k c 



where < £ h < By its definition, ^ h represents the amount of electricity that can be supplied in 
the market from coal generators whose heat rate is smaller than or equal to h c . 

Proposition 3. Given I — {c,g}, if, under Q, the random variables log {S C T ) and log(Sf.) are jointly 
Gaussian distributed with means \\>c and jig, variances (J^ and (X? and correlation p, then, for 



-i 

Explicit formulae for cases of h c outside of this range are also available, but not included here. 
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t € [0, T], the price of a dark spread option with maturity T is given by 
(14) V t =exp{-r(T-t))lj o v w , 2 (D,F t ) (j> d (D) dD 

f<- ,mas(f*,5*) 

,. . viow i (D,F t ) tj>d(D) dD + I v mid3 (D,F t )<j) d (D)dD 

«nin(f*+$*f e ) /•!''+ 

+ / ,„ ts v W .2,, + (£>,F,)fe(D)dD+ / v mWi i(D,F,)^(i>)dD 

mK(| e ,f'+5*) /•! ) 

v h igh,2 (D,F t ) <j) d (D) dD + vu g h,\ (D,F t ) <j) d {D) dD } , 



where the integrands are given in Appendix\B\and discussed in some detail. 

Proof. As in the proof of Proposition^ by iterated conditioning, for t e [0, T], the price of the dark 
spread V, is given by 

V, : = exp (-r(T - 1)) E® [ (P T - h c S c T ) + \ ft] 

= exp(-r(7'-r))E Q \E Q [(b(D T ,S T ) - h c S' T ) + \^W ft w ] ft 

Again we write the outer expectation as the sum of integrals corresponding to the different forms 
the payoff can take, since the functional form of b is different for Dj lying in the intervals [0,£'~], 
[<^'-,f'+], and [^ ,+ ,^]. In addition the functional form of the payoff now depends on whether 
Dj < E, h or Dj > £, h and on the magnitude of t, h relative to f c and t, 8 . Therefore, the first case is 
subdivided into the intervals [0,min(| s ',^' ! )], and [min(f g ,^ A ),f *-]; the second case is subdivided 
into [fS max(f [max(^,^),min(|« + §M c )], and [min(|* + ^*,| c ),|''+]; the third case 
is subdivided into [| i +,niax(| c ,|* + ^ A )], and [max(f^|* + |*),|]. 

The integrands v in ( fl4] > are obtained by calculating the inner expectation for each demand inter- 
val listed above, in a similar fashion as in Proposition [2] □ 

Note that ( fl4| i requires seven terms in order to cover all possible values of h c within the range 
given by ( |13) , as well as the two cases c = i + and c = z'_ . However, only five of the seven terms 
appear at once, with only the second or third appearing (depending on h c ^ exp(k c +m c ^ 8 )) and only 
the fifth or sixth (depending on h c ^ exp(k c + m c (^ c — B, s ))). These conditions can equivalently be 
written as % h | s and § A «g | c - f Notice that if c = /_, we deduce that £, h <l s and > | c - | g 
irrespective of while for c = z' + several cases are possible. 

Similar to the analysis of the forward contract earlier, if demand is assumed to be deterministic, 
then the spread option price is given explicitly by choosing the appropriate integrand from Proposi- 
tion [2] To now obtain a convenient closed-form result for unknown demand, we extend our earlier 
notational tool for combining Gaussian distribution functions. For any integer n, let 



Xu X\2 ■■■ X\ n 
X 2 \ X 2 2 ■■■ X 2n 



,y\P ) =E[* 2 ( Xl ''^ ; ' C) )- 4> 2(^2i,3';p)]- 

;=1 



In addition, we introduce the following notation to capture all the relevant limits of integration. 
Define the vector a := {a\,... ,ag) by 

a:=~[(0, §*A§* I''-, |«v§* § c A + f<'+, | C V(|« + ^), |)-tfi]. 

Of/ 
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Notice that the components of a are in increasing order and correspond to the limits of integration 
in equation ( fT4| ). In the case that c = i + , all of these values are needed, while the case c = z"_ is 
somewhat simpler because 03 = 04 and 05 = ag (since by ( fT3| l, E, h < f c ). However, the result below 
is valid in both cases as various terms simply drop out in the latter case. 

Corollary 4. In addition to the assumptions in Proposition^let demand at maturity satisfy 

D T = max (0,min (%,Xt)) , 

where Xj ~ N(jX d , C7?) is independent of St- Then for t G [0, T], the price of a dark spread is given 
explicitly by 

R c {H d ,Q)-m 2 c aj m c a d 



V t = exp(-r(T-t)) lb c (n d ,F t c )exp 
+ b c (li d -Z 8 ,F t c )exp 

+ ^(Md-| e ,F/)exp 



(t aJ )^ x2 ( 




a 3 




04 


«2 





fl 8 


«6 




a-/ 


fl 5 



# c (ft* ~ I s , l g ) + m z c aj -m c a d \ 



2x1 



<?cJ 



Gc,d 

-m 2 g aj -m g a d ^ 



-h c F t c <t> 2 * 3 



a-j 05 03 

flg fl4 fl2 



«8 

R c ((loghc - j3 - YUa) /a g ) _ -ycT rf 



g u g,r 



+fc C j(Md,F,)exp(T}) <^ 4> 



.2x2 



f C A3 
«4 «2 



-i? c (na ,0)-a g a 2 + ym c aj -m c a d 



4> 



2x2 



flg (26 
07 05 



4> 



2x3 



+ 4) 2xl 



fl7 05 03 
flg fl4 fl2 



«8 



Gc,d Gc,d 

Oc,</ ' Oc,d y 

/?g (tfd ~ I c , I c ) + « c g 2 - ym g <J 2 m g a d \ 
a g4 ' <*g,d J 

-r c ((log// - p - r MdVgg) - « g g 2 - y 2 g|/«g . ygrf 

a„ y ' a,, (7, 



*l(-«8)E*i(l'»^)*l 

(6/ 



/i CJ F, c (1 - *j (fl 7 ) + *! (fl 6 ) - *! (as)) 1 , 



where 



^(z):=z + log(F/)-log(/;')-icJ 2 flnfl 7 a 2 , := y^/a, 2 + a 2 



Proof. All terms in ( fi4] i have the same form as those in Proposition|2]for forwards: demand appears 
linearly inside each Gaussian distribution function and in the exponential function multiplying it. 
Hence, applying Lemma[T]and simplifying lead to the result in the corollary. □ 



Remark 7. Under the extended model introduced in f 3.3 to capture spikes, spread prices are given 
by the same expression as in Proposition|4]plus the following simple term^] 



exp ( m s (ji d - 1) + -m 2 s a 2 d ) 4>, 



m s a d 



*1 



o d 



^We require only two of the four extra terms in Remark^due to our assumption on h c in {13) , which guarantees that for 
the spike regime, the option is always in the money, while for the negative price regime, it never is. Hence, if we were to 
consider put spread options instead of calls, the other two terms would be needed instead. 
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Bid curves 


Fuel price processes 


Demand 


Rate 


ki nii ^ ' 


Ki Vi Aj s' 




r 


2 1 0.5 


1 0.5 log(10) 10 


0.5 0.2 






Table 2. Parameters used throughout §6 (for ; e I) 



6. Numerical Analysis of Spread Option Prices 

^on 



In this section, we investigate the implications of the two-fuel exponential stack model of \ 3.2 
spread option prices and power plant valuation, as compared to two common alternative approaches. 
We analyse prices for various parameter choices and option characteristics, as well as fuel forward 
curve scenarios. 

Recall that the closed-form spread option prices given by Corollary [4] required no specification 
of a stochastic model for fuel prices, but instead imposed only a lognormality condition at maturity 
T. However, for the purpose of comparing prices across maturities and across modeling approaches, 
we select a simple example of fuel price dynamics consistent with Corollary|4] We assume coal (S,) 
and gas (Sf ) follow correlated exponential OU processes under the measure Q. i.e., for i £ I and 

te[o,T], 

(15) d(logSj) = K,(X I -(logS i t ))dt + v i dW t i , 5(, = y 

where d (W c , W g ) t = pdt. As gas and coal are treated identically in the bid stack model, we consider 
the symmetric case in which all coal and gas parameters are equal j^] including both the parameters 
in ( fT5j ) and in the exponential fuel bid curves, defined in |5]). All are listed in Table [2] with the 
exception of p, which we vary throughout our analysis. All prices are calculated for time t = 0. 
Note that for a given maturity T, the parameters jx c , fl g , O c , and p appearing in f|4]and f|5]are 
related to those in ( fT5| ) by the following standard results (for i 6 /): 

IU = i'exp(-K-,T)+A 1 (l-exp(-K- i r)), 
of = ^(1- exp (-2*7-)), 
P G ' a s = ^(l^xpKr-)C s r)). 

K c -f - Kg 

As for fuel prices, recall that no particular process is required for electricity demand in our model. 
Typically driven by temperature, demand is often modelled as rapidly mean-reverting to a seasonal 
level. Hence, in our examples, we assume (X t ), as introduced in Corollary|4j to be an independent 
OU process, with a high value for mean-reversion speed (e.g. 100 or more). For our aim of pricing 
options with maturities of several months or even years, the values of these parameters are insignifi- 
cant, as Dj is always well approximated by its stationary distribution. Hence, in Table|2]we list only 
the values and a ( /- We also assume the interest rate r = throughout. This completes the base 
parameter set to be used throughout this section unless otherwise stated. 

Next, we introduce three scenarios designed to assess the role of observed fuel forward curves 
and corresponding implications for bid stack structure. Fuel forward curves reveal crucial informa- 
tion about the probability of future merit order changes. For example, if bids from coal and gas 
are currently at similar levels (as in our parameter set) but one fuel is in backwardation (forwards 



"'While it is of course not realistic for gas and coal to have identical prices (typically differ by a magnitude of about ten 
and trade in different units), it is plausible that after adjusting for different heat rates, the coal and gas bids could indeed 
coincide and have similar volatility, make our symmetric case less hypothetical. 
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decreasing in T), while the other is in contango (increasing in T), then the future dynamics of power 
prices (under Q) should reflect the high chance of the coal and gas bids separating. We compare the 
following scenarios: 

(I) No fuel forward inputs (forward prices implied from ( fT5| ) 
(II) Gas in contango; coal in backwardation (linear with stepsize of 0.2 per month) 
(III) No fuel forwards, but gas bids above coal (X c = log(s c ) = log(7), X g = log(i 4 *) = log(13)) 
Note that in Scenario II, fuel forward curves are inputs assumed to be observed from the market, in- 
stead of being generated by a model. The standard approach to resolving the inconsistency between 
the market and the model in ( fTB) is to calibrate to fuel forwards for each T via a shift in the mean 
level jj,, (or more formally via a time dependent long-term mean A,). 

6.1. Spread Option Price Comparison. To test our model's prices for spark and dark spread op- 
tions, we compare with two other typical approaches to spread option pricing: Margrabe's formula 
(cf. EBI ) and a simple cointegration model (cf. ifTTI . (15 1 for discussions of cointegration between 
electricity and fuels). 

6.1.1. Margrabe's Formula. Assume that under the measure Q, the electricity price Pj and fuel 
price S l T , i £ I, are jointly lognormal, with correlation p p ,-. Writing fi p and a p for the mean and 
variance of log(Pj), then for t <E [0, T\, the price of a spread option with payoff ( p~2] > is given by 

T „ (\og(F t p /hiFi)+al J2\ f log (Ff/hiFf) - (7,2/2 V 
Vr=exp(-r(r-0) Ff*! f * V 'J [ '^—\ -hfi*i f &K ' J 1 , 

where a p l ■ — <J p — 2p ; , ,a p C7/ + <7, 2 , and all other notation is as before. 

6.1.2. Cointegration Model. Let Yj be an independent Gaussian random variable under Q, with 
mean jj, y and variance a y . Then for constant weights w c ,w g > (the cointegrating vector), we 
define Pj by 

P T := w c S c T + w g S k j+Y T . 
No closed form results are available for spread options, so prices are determined by simulation. 

6. 1.3. Comments on Comparison Methodology. In order to achieve a sensible comparison between 
the stack model and either Margrabe or the cointegration model, the mean and variance of Pj should 
be chosen appropriately, for each maturity. For a single fixed T this simply requires choosing pa- 
rameters jip and Op (in the case of Margrabe) or fly and a y (in the case of the cointegration model) 
to exactly match the mean and the variance of Pj produced by the stack model. As we shall often 
compare models simultaneously across many maturities, we use a variation of this idea, finding the 
best fit of an OU process for log(P f ) or (Y t ). For Margrabe, 

d(logP f ) - K p {X p - (logP,)) df + V p dW t p , P =p , 

where d(W p ,W) = Ppjdt for i e /, while for the cointegration case (Y t ) is an independent OU 
process. In all cases we utilize the closed form expression for the variance of Pj in the stack model, 
given in Appendix |A| 

Finally, to reflect the symmetry between coal and gas in our parameter set, in the cointegration 
model we set w c = w g = 1 /2exp(£ c . +m c (/x</4 ' )), such that the power price has equal dependence on 
each fuel, with a price level linked to the most likely marginal bid levels of coal and gas. Note that 
since both underlying fuels appear in the cointegration model like in the stack model, a comparison 
across values of fuel price correlation p is quite natural. On the other hand, the correlation parameter 
which we vary in the Margrabe approach is p p which correlates electricity and fuel. Hence the 
choices of p or p Pj , in the following plots are not perfect comparisons, but rather illustrate the role 
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of correlation in generating a range of prices for each model. Due to variations across maturities, 
there is no direct link between a chosen value p and an appropriate p p j short of estimating these 
from data. 
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Figure 6. Spread option prices against ly for different correlations and demand levels 



6.2. Spread Option Parameter Analysis. 

6.2.1. Spread prices versus heat rate hi. For a fixed maturity T = 1, we plot dark spread option 
prices as a function of heat rate h c (over a range corresponding to (fT3)>). Figure 6a illustrates the 



Scenario I results for each of the three models considered, for two different correlation parameters, 
p = ±0.8. In all three models, negative correlation logically raises the option price, as it increases 
the volatility of the spread. The stack model generally predicts lower spread option prices than Mar- 
grabe, and a smaller gap between different correlation levels. This is because the strong structural 
link keeps long-term levels of power and gas close together, thus narrowing the spread distribution 
relative to the weaker case of correlated Brownian Motions. The gap between Margrabe and the 
stack model widens with negative correlation, while for extreme positive correlation (not plotted), 
Margrabe can underprice the stack model. Another way of understanding this phenomenon is to no- 
tice that the bid stack automatically imposes a positive dependence structure between electricity and 
its underlying fuels, which can only be somewhat weakened or strengthened by varying fuel price 
correlation through p. The cointegration approach shares this characteristic, and therefore prices 
much closer to the stack model than Margrabe, but still somewhat higher for p = 0.8. 



Unlike for Margrabe, Figure 6b reveals interestingly that the dependence on correlation p does 
not hold strictly for the stack model, as exceptions can be found. Considering a market with very 
low demand (ji c i = 0.3, G c { — 0. 12), we observe that the relationship with p is reversed for low values 
of h c . Since demand is so low, the price is typically set by the cheaper fuel and hence only in the 
money for coal bids generally below gas. Hence for negative p, we do not receive the typical high 
payoffs from high gas and low coal states, and get zero payoff when coal moves above gas in the 



stack. Figure 6b also confirms that this reversal is most pronounced if coal bids are shifted to be 



below gas bids, as described by Scenario III. 

6.2.2. Spread prices versus maturity T. We next investigate spread option prices against maturity T, 
and for this purpose fix h\ — e\p(k, ■ + ;«;<i;'/2) in option payoffs, matching the median heat rate in the 
market. We again compare several correlation levels for all three models, and now include Scenario 



II to test the impact of fuel forward curves. In Figure 7a (Scenario I), spread options are typically 



increasing in maturity as expected, and flatten out as the price processes approach their stationary 



distributions, with an ordering of the three models resembling Figure 6a More interestingly, in 
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Figure 7. Option prices against T for different correlations and fuel forward scenarios 



Scenario II (Figure 7b I, longer term spark spread options drop significantly in value in stack model, 
are thus greatly overpriced by Margrabe, and significantly overpriced by the cointegration model as 
well. As the gas forward curve is now in contango, while coal is in backwardation, coal will almost 
always be below gas in the future bid stack, especially for very long maturities. Hence, a spark 
spread option has relatively little chance of being in the money, as this would require unusually high 
demand. This is a good example of a dependency which cannot be captured by Margrabe or other 
reduced-form models, but is automatically captured by the merit order built into the stack model. 
Moreover, fuel forward prices are direct inputs into our expressions for spread options, avoiding the 
need for an additional calibration step to first match observed fuel forwards, as is the case for the 
other approaches. 
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Figure 8. Analysis of impact of matching mean and variance of Pj for Scenario II 
So far all plots have assumed that both the mean and variance of the power price distributions 



are matched in all three models, via the procedure described in Section 6.1.3 One might question 
whether this is realistic. In practice, we only have history (and possibly observed forward curves) 
to calibrate each model, and thus should not be borrowing extra information about the future from 
the stack model's structure when calibrating the other approaches. While matching the mean is 
reasonable as it is analogous to matching observed power forwards, matching the variance is less 



justifiable. In Figure 8a we compare spark spread option prices for Margrabe and the cointegration 
model in Scenario II (as in Figure 7b but now p = 0) for three different calibration assumptions: 
full matching as earlier; matching means but not variances; matching neither means nor variances. 
Here 'not matched' implies that means and/or variances are instead fitted to Scenario I levels (a 
proxy for history). We note significant differences between all cases. Failing to match the mean 
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implies a greater overpricing of spark spreads in this case, while failing to match the variances acts 
in the opposite direction here, lowering the price since the forward looking variance (implied by the 
stack in Scenario II) is higher than the variance in Scenario I (see Figure |8b]>. While other scenarios 
could lead to different patterns, it is clear that significant price differences can occur due to the 
likely changes in the merit order. In Margrabe, no information is transmitted from fuel forward 
prices to the distribution of Pj, while in the cointegration model limited information is transmitted, 
since the relative dependence on coal and gas is fixed initially by w c ,w g , instead of dynamically 
adapting to fuel price movements (and demand). In contrast, the stack model produces highly state- 
dependent power price volatility and correlations reflecting known information about the future 
market structure. 
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(C) p"" p surface (Scenario I) (D) p' mp surface (Scenario II) 

FIGURE 9. Implied correlation analysis for various parameters and scenarios 



6.3. Implied Correlation Analysis. We next analyse 'implied correlation' p ' , meaning the value 
of p p j for which Margrabe's formula reproduces the stack model price. As Figures|6]and|7]suggest, 
for high (positive) values of p in the stack, it may be impossible for Margrabe to reproduce the price, 
for any p p .[ £ [—1,1]. In such cases, implied correlation does not exist. However, p™ p typically 
exists for most values of p, and can be understood as a convenient way of measuring (or quoting) 
the gap between Margrabe and the stack model price. 

In Figures 9apb we investigate implied correlation 'smiles' (against h c ) for a dark spread option 
in Scenario I and with p = 0. In Figure 9a we first vary the relative capacities of coal and gas (with 
^ = 1 throughout). In all cases Margrabe overprices the spread since p' p m f > 0, but the difference 
is much larger when coal is the dominant technology. As we approach the case of a dark spread 
in a fully gas driven market (f g near 1), Margrabe approaches the stack price (p imp near 0). In 
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Figure 9b we assume ^ c = ^ g = 0.5, but instead vary jj, c i (with G c \ now 0.12). We see that the 



implied correlation has a slight downward (upward) skew if demand is high (low), and a fairly 
symmetric 'frown' for jij = 0.5. Figures 9cpd plot implied correlation as a function of both h and 
T for Scenarios I and II. When given fuel forward curves as inputs (Scenario II), we can observe a 
distinctive tilt in the implied correlation surface for long maturities. 

6.4. Power Plant Valuation. We conclude this analysis with an investigation into the bid stack 
model's predictions for power plant valuation under our chosen scenarios. A generating unit of fuel 
type i E I, with heat rate hj, can be approximated as a sum of spread options on spot power (cf. 
03D))- Letting {Tj : j E J} represent all future hours of production over the plant's life, the plant 
value (VP t ) for t E [0, T] is then 



VP, = £exp (-#■(?)■- f))E< 
jeJ 



Prj-hiSir. 



While this approximation technique ignores complicated operational constraints, it is consistent with 
our approximation of the electricity price setting mechanism itself, since a plant bidding at cost every 
day receives exactly this string of payoffs in our model. 
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Figure 10. Power Plant Value (3yrs, 1000MW) as a function of fid for various scenarios 



In Figure 10 we value a 1000 MW power plant with a life of three years in the Margrabe case 
and stack model, exploiting our closed-form formulas for reasonable computation time. Instead 
of making an arbitrary assumption about the important periodicities of power demand (which vary 
from market to market) we use a fixed mean fid for all hours in each calculation, but investigate 
the resulting plant value as a function of We also fix hi = exp(fe,- + m,f */2) throughout. In 
Scenario I (Figure |10a|i, Margrabe prices consistently higher than the stack model as expected. 



Figure 10b considers Scenario II where fuel forward prices pull future gas and coal bids in opposite 
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directions, and Margrabe matches the distribution of Pj only from history. Remarkably, for long 
enough maturities and high enough jj, c /, Margrabe sometimes underprices a coal power plant relative 
to the stack model. Here the stack captures that gas is likely to be the marginal fuel in the future, 
with coal plants operating near full capacity. 

Another interesting case to consider is Scenario III, in which there is very little overlap between 



coal and gas bids to begin with, and coal is likely to remain below gas in the merit order. Figures 10c 



lOd reveal the result of this change. Unsurprisingly, for a gas power plant (Figure lOd i the deviation 
between the stack model and Margrabe is large, since the gas plant has little chance of being called 
upon to produce power. On the other hand, the difference between Margrabe and the stack model 
is much less for the coal plant, and the stack model price appears to converge to Margrabe for high 
demand (similarly to Figure 9a for high £, g ). The reason for this is that when coal is always below 
gas in the stack and demand always high, then the power price can be approximated by the gas stack 
alone. Hence, power price should be close to lognormal and the correlation between power and coal 
close to that of gas and coal. Furthermore, as there is no mismatch of mean or variance (as there 
was in |10b| i, under such a scenario, Margrabe's formula should give a very similar price to the stack 
model. 

Finally, it is important to remember the heavy-tailed nature of most electricity spot prices, which 
makes a lognormal distribution for Pj highly questionable. As discussed in Section 3.3 the bid 



stack model allows for a straightforward extension to capture spikes (or negative prices) consistently 
without limiting the availability of derivative pricing formulae. To illustrate the impact of spikes, 
Figur eJIT] plots the same scenarios as in Figures 10a and lOd except using our extended stack 
moderfwith parameter m s = 50 (and a c i reduced from 0.2 to 0.1). The impact of spikes is to shift 
the stack model valuation closer to Margrabe, but only in the case of high /i c /. This is intuitive of 
course because spikes are only present for high enough demand. Since the stack model produces a 
heavier-tailed distribution for Pj than the lognormal, very high ji c i can lead to stack model valuation 
even above Margrabe. 
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Figure 1 1 . Power plant value with stack model extended to include spike regime 



7. Conclusion 

The valuation and hedging of both physical and financial assets in electricity markets is a com- 
plex and highly state-dependent challenge, particularly for medium to long term problems in markets 

^This implies a spike of about $150 in the event X t — 1.1 f, and so is a fairly conservative choice. Note also that we still 
match the mean of the Margrabe formula (i.e., calibrate to the stack model's electricity forward curve), but only match the 
variance to the case of no spikes, since otherwise the fit to the new variance (with spikes) becomes infeasible for large fi^. 
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driven by multiple underlying fuels. As we have seen in the examples of Sj6j it is important to be 
aware of the current merit order and resulting bid stack shape, as well as likely or possible changes 
to this order in the future. Observed fuel forward curves can provide valuable information for this 
purpose, but cannot be incorporated easily into traditional reduced-form models for power prices. 
On the other hand, a structural approach maintains a close link with the physical characteristics of 
the electricity market, allowing for the inclusion of a variety of forward looking information, such 
as demand forecasts, or changes in the generation mix of the market, a pertinent issue in many coun- 
tries nowadays. The piecewise exponential bid stack model proposed here achieves this link, while 
crucially retaining closed-form expressions for forwards and spread options, as presented in Sec- 
tions [4] and [5] In this way, it enjoys the benefits of a simple reduced-form model, while mimicking 
the complex dependence structure produced by a full production cost optimization model, for which 
derivative pricing is typically a computationally infeasible task. Furthermore, the availability of ex- 
plicit expressions for forwards and options is highly beneficial for many other practical purposes, 
such as calibration to observed market quotes, the calculation of Greeks and for risk management 
applications (e.g. simulating price changes for a portfolio of physical assets). Although we focused 
primarily here on a simple coal and gas based market, we believe that our general framework has 
the flexibility to be adapted to various market conditions, as illustrated for example by our simple 
extension to capture extreme spikes at times of high demand, an important feature of most power 
markets. As we have shown through many examples, the delicate interplay between demand, capac- 
ity, and multiple fuel prices is embedded into our approach, thus providing an intuitive framework 
for understanding complicated interdependencies, while also helping to bridge the prevalent gap 
between mathematical tractability and the economics of supply and demand. 



Appendix A. Moments and Covariances 
If demand at maturity satisfies 

D T = max (0,min (|,X r )) , 
with Xj ~ N(jld,Od) independent of S7-, then for t g [0, T], the n-th moment of Pj is given by 
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where j = I\{i}, Si = (— 1) {*-'+} as before and 

Rf\^j) :=k j +m£ j -k i -m&+\og{F t j )-\og{F t i ) - {n- l -)af +np<no h 

ri (n) --=Y - <W 2 ) , and ^ := \ („ 2 - n)of + ^mfal 

A more general formula for E[Pj (S C T )" c (S^)" s \ ,^ t ] can be obtained similarly, allowing us to calculate 
for example covariances between electricity and fuels. 

Appendix B. Terms in Spread Formula 

The integrands in Proposition|4]correspond to the dark spread price in various demand intervals, 
in the case of a given or known demand value. These terms resemble those for forwards in Proposi- 
tion [2] but can be categorized by whether the option is always, never or sometimes in the money for 
each case. For the last of these three cases, we have 

y lowJ (§,*)= b c (<^ £ ')4>i (* c (S,0)/a) -h c x c ^ {R c /a) 

+b cg exp (a^) [1 - <£, (R c ($ , 0)/a - a ? a) - *j (tf, /a - a c a)] , 

VMgh.2 (| ,x) = - f (-J? g (4 - | c , | c ) /a) - ^4>i (* e /a) 

+ ^ (4 ,x) exp (- CC £^1\ ( Rg (| _ ^ | C ) /a _ acty ) _ 4 , 1 (* g j a _ a cty )] , 

and v midAc - = viow.i, Vmid,2, g = v high , 2 , where = (logh c - j8 - Y^)/ a g- Then 

Vmid.l =/mid(4> x )- ft c^ C , Vhigh.l (|,x) = /high (4 : x ) - 

correspond to cases that are always in the money, while the following are never in the money: 

Vlow,2 (4: x ) = Vmid,3(4> x ) =0 
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